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Abstract. The classical vacuum gap model of Ruderman & Sutherland, in which 
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spark-associated subbeams of subpulse emission circulate around the magnetic 

, axis due to the E X B drift of spark plasma filaments, provides a natural and 

plausible physical mechanism of the subpulse drift phenomenon. Moreover, this is 

q . the only model with quantitative predictions that can be compared with observa- 

5h ' 

■ tions. Recent progress in the analysis of drifting subpulses in pulsars has provided 

C/3 . 

a strong support to this model by revealing a number of subbeams circulating 

around the magnetic axis in a manner compatible with theoretical predictions. 
• *~ i ■ 

However, a more detailed analysis revealed that the circulation speed in a pure 

\—( ' 

vacuum gap is too high when compared with observations. Moreover, some pul- 
sars demonstrate significant time variations of the drift rate, including a change of 
the apparent drift direction, which is obviously inconsistent with the E x B drift 
scenario in a pure vacuum gap. We attempted to resolve these discrepancies by 
considering a partial flow of iron ions from the positively charged polar cap, co- 
existing with the production of outflowing electron-positron plasmas. The model 
of such charge-depleted acceleration region is highly sensitive to both the critical 
ion temperature T, ~ 10 6 K (above which ions flow freely with the corotational 
charge density) and the actual surface temperature T s of the polar cap, heated 
by the bombardment of ultra-relativistic charged particles. By fitting the obser- 
vationally deduced drift-rates to the theoretical values, we managed to estimate 
polar cap surface temperatures in a number of pulsars. The estimated surface 
temperatures T 3 correspond to a small charge depletion of the order of a few 
percent of the Goldreich- Julian corotational charge density. Nevertheless, the re- 
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maining acceleration potential drop is high enough to discharge through a system 
of sparks, cycling on and off on a natural time-scales described by the Ruderman 
& Sutherland model. We also argue that if the thermionic electron outflow from 
the surface of a negatively charged polar cap is slightly below the Goldreich- Julian 
density, then the resulting small charge depletion will have similar consequences 
as in the case of the ions outflow. We thus believe that the sparking discharge 
of a partially shielded acceleration potential drop occurs in all pulsars, with both 
positively ("pulsars") and negatively ("anti-pulsars") charged polar caps. 

Key words, pulsars: general: plasmas - pulsars: individual: PSRs B0943+10, 
B0809+74, B0826-34, B2303+30, B2319+60, B0031-07 

1. Introduction 

The phenomenon of drifting subpulses has been widely regarded as a powerful diagnostic 
tool for the investigations of mechanisms of pulsar radio emission. The drifting subpulses 
change phase from one pulse to another in a very organized manner, to the extent that 
they form apparent driftbands of duration from several to a few tenths of conscqutive 
pulses. Typically, more than one drift band appears, and the separation P3 between them 
measured in pulsar periods Pi ranges from about 1 to about 15 (Backer 1973. Rankin 
1986 this paper). There are typically two or three approximately equidistant subpulses 
in each single pulse, separated by P2 degrees of longitude. Thus, the observed drift rate 
Do = P2/P3 degrees of longitude per pulse period Pi , provided that P3 is alias- free or 
alias-corrected (e.g. van Leeuven 2003; see also Gil & Sendyk 2003 and references therein 
for review). The subpulse intensity is systematically modulated along drift bands, either 
decreasing (typically) or increasing (seldom) towards the edge of the pulse window. In 
some pulsars, however, only periodic intensity modulations are observed, without any 
systematic phase change. These pulsars were identified as those in which the line-of- 
sight cuts through a beam centrally (Backer fl973ll . thus showing a steep gradient of the 
polarization angle curve (e.g. Lyne & Manchester I1988|l . On the other hand, the clear 
subpulse driftbands are typically found in pulsars associated with the line-of-sight grazing 
the beam, thus showing relatively fiat position angle curve (Backer [19731 Rankin H 9861 
Lyne & Manchester I1988f) . The observed periodicities related to patterns of drifting 
subpulses are independent of radio frequency, thus excluding all frequency dependent 
plasma effects as a plausible source of the drifting subpulse phenomena. 

The observational characteristics of drifting subpulses described briefly above suggest 
unequivocally the interpretation of this phenomenon as a number of isolated subbeams 
of radio emission, spaced more or less uniformly in the magnetic azimuth, and rotating 
slowly around the magnetic axis. The most spectacular confirmation of this interpretation 
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was recently presented by Deshpande & Rankin (1999, 2001; DR99 and DR01 hereafter) 
and Asgekar & Deshpande l|2001|l . who performed a sophisticated fluctuation spectra 
analysis of single pulse data from PSR 0943+10 and detected clear spectral features 
corresponding to the rotational behaviour of subpulse beams. If each subbeam completes 
one full rotation around the magnetic axis in P3 pulsar periods Pi, then N = P3/P3 
is the number of circulating subbeams. The primary periodicity P3 is relatively easy to 
measure, either by eye-inspection or by finding a corresponding frequency fa = I/P3 
in the intensity modulation spectrum (although in some cases it requires an aliasing 
resolving - see Gil & Sendyk 2003, van Leeuwen et al. 2003, DR01, DR99). However, 
measuring or estimating the circulational (tertiary) periodicity P3 is much more difficult, 
since detecting a corresponding feature in the fluctuation spectrum at low frequency 
/3 = I/P3 requires an extraordinary stability of intensity patterns over a relatively long 
period of time. So far, such a low frequency feature was found only in the fluctuation 
spectrum of PSR 0943+10, Asgekar & Deshpande (2001; see their Figs. 1 and 2) detected 
clear peak at fs — 0.027/Pi, corresponding to circulational period P3 s» 37Pi (see Gil & 
Sendyk 2003, for more detailed discussion). Moreover, DR99 & DR01 detected sideband 
features near the high frequency feature /3 = I/P3 (separated from it by 0.027/Pi), also 
clearly associated with the rotational cycle P3 ~ 37Pi in PSR 0943+10. Van Leeuwen et 
al. (2003) analyzed drifting subpulses in the well known pulsar PSR 0809+74 and argued 
that P3 > 150Pi and N > 14 in this pulsar. Recently, Gupta et al. I|2003|l analyzed 
drifting subpulses in PSR 0826-34 and found that P 3 w P 1; P 3 = 14P[ and N = 14 in 
this pulsar (see Sect. 4.3 for some details of this analysis). We will use the observationally 
deduced values of circulational periodicities P3 of these three pulsars (and a few others) 
later on in this paper in order to estimate basic parameters of the polar cap physics. 

The frequency independence of the drifting periodicities P3 and P3, as well as sim- 
ilarities of the drifting subpulse patterns at different radio frequencies strongly suggest 
that the radiation subbeams in the emission region reflect some kind of a "seeding" 
phenomenon at or very near the surface of the polar cap (rather than some kind of mag- 
netospheric plasma waves; e.g. Kazbegi et al. 1996). As DR99 emphasize, the results of 
their analysis of drifting subpulses in PSR 0943+10 appear fully compatible with the 
Ruderman & Sutherland 1)19751 RS hereafter) E x B drift model, although their analysis 
is completely independent of this model. In this paper we provide further support to this 
natural model, whose original version we review shortly in Appendix A. 2. In Sect. 2 we 
discuss a modified version of RS model, allowing a partial ion or electron flow from the 
polar cap surface, coexisting with the magnetic electron-positron pair plasma produc- 
tion. In Sect. 3 we discuss a thermostatic regulation of the polar cap and estimate actual 
surface temperatures. In Sect. 4 we calculate the predicted circulational periodicities and 
compare them with the observationally inferred values for a number of pulsars. Finally, 
we give a summary of our results in Sect. 6. 
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The original version of RS model can be applied only to pulsars with a positively 
charged polar cap ("pulsars" in RS terminology), i.e. with n • B < 0, which is usually 
considered as a deficiency of the model. We propose a natural solution for the other 
"half" of neutron stars ( "antipulsars" in RS terminology), with fl ■ B > 0, in Sect. 2.2. 
RS assumed a strong binding of iron ions, which therefore could not be released from the 
polar cap surface by thermionic and/or field emission. In this paper we argue that even if 
the iron ions (or electrons) are marginally bound within the surface, then the centrifugal 
outflow of charges through the light cylinder results in the creation of an acceleration 
region just above the polar cap surface. The residual potential drop is strong enough 
to be discharged by the magnetic creation of electron-positron (e _ e + ) pairs that form a 
system of isolated plasma filaments (sparks) , which in turn produce a system of isolated 
plasma streams flowing along dipolar magnetic field lines and radiating spark-associated 
coherent subpulse radio emission at higher altitudes (Kijak & Gil 1998; Melikidze et al. 
2000). 

The important feature of the RS model is an inevitable E x B drift, which makes 
the spark plasma filaments rotate about the symmetry axis of the surface magnetic field. 
Gil et al. (1993) argued that this circumferential motion of sparks is manifested by conal 
structure of pulsar beams (Rankin 1983, 1986). RS adopted a pure axial symmetry of 
a star-centered global dipolar field, although they implicitly assumed significantly non- 
dipolar radii of curvature TZq ~ 1 (see Appendix A.l.) required by conditions of the 
magnetic pair production. Both the spark characteristic dimension, as well as the distance 
between adjacent sparks should be about the height ft, of a quasi-steady vacuum gap (RS, 
Gil & Sendvk EZoTJUl hereafter GS00). The speed of the E x B drift motion around the pole 
is Vd — cAE±/B s cm s — 1 , where AE± is the component of the electric field caused by the 
charge depletion in the acceleration region. A prominent subpulse drift can be observed 
when the line-of-sight grazes the pulsar beam, which corresponds to peripheral sparks 
drifting at a distance d ~ r p — h from the pole (see Appendix A. 2. for details). Each spark 
completes one full rotation around the magnetic axis in P3 = 2ird/vd seconds (called the 
tertiary periodicity). According to Eqs. (A. 3) and (A. 4), for d ~ r p — h and 77 = 1 (pure 
RS vacuum gap) the tertiary periodicity P3/P1 = [(r p /h) — 1] and the azimuthal drift 
rate D r = 360°/p3 = 360°/[(r p //i) — 1]. It is important to emphasize that the value of 
D r can be deduced from observations of drifting subpulses only if the value of P3 can be 
measured/estimated. On the other hand, D r can be theoretically estimated if the value 
of the ratio r p /h is known. In the case of PSR 0943+10, the value of r p /h ~ 7 (see GS00 
and Gil et al. 2002b; GMM02b hereafter;) and P 3 = 37.35 Pi. Thus, the vacuum drift 
(RS) periodicity P3 ~ 6 Pi is about six times shorter than the observed value P3 s» 37Pi, 
and, consequently, the drift rate D r = 360° /P3 is about 6 times too high (see also Gil & 
Sendyk 2003). In PSR 0809+74, which is another pulsar for which P3 could be estimated, 
this discrepancy is much larger. In fact, van Leeuwen et al. 2002 demonstrated that the 
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observationally deduced value of P3 exceeds 150Pi, while the RS value of P3 is about 
5Pi. Similarly, in PSR 0826-34 the RS model gives P 3 = 5Pi, while Gupta et al. (2003) 
deduced from the analysis of drifting subpulses that P3 ~ 14Pl in this pulsar. 

Therefore, in pulsars for which the azimuthal (intrinsic) drift rates D r can be mea- 
sured/estimated (see Table 1), they turn out to be a few to several times lower than 
those predicted from RS model. In other words, the pure vacuum gap drift is too fast as 
compared with observations. Moreover, in some pulsars a time variable drift rate is ob- 
served, including reversals of apparent drift direction in a few cases. These observational 
features are inconsistent with the RS model, which otherwise provides a quite natural 
and plausible physical mechanism of the subpulsc drift phenomenon (not to mention 
that this is the only quantitative model that can be compared with observations). In this 
paper we attempt to resolve these discrepancies within a more general model of the inner 
acceleration region, involving a partial flow of iron ions (or electrons) due to the thermal 
emission from the polar cap surface, heated to high temperatures by sparking discharges. 
Such generalization of the pure vacuum gap model of RS was first proposed by Cheng & 
Ruderman (1980, CR80 hereafter; see also Usov & Melrose 1995, 1996). However, CR80 
suggested that even with ions included in the flow, the conditions above the polar cap are 
close to a pure vacuum gap. We, on the contrary, argue in this paper on both theoretical 
and observational grounds that the quasi-stationary discharge conditions can be estab- 
lished, even if the ion (electron) flow exceeds 95% of the Goldreich & Julian (1969; GJ 
hereafter) charge density. Nevertheless, the remaining acceleration potential drop is high 
enough to discharge through a system of sparks, as originally proposed by RS. The im- 
portant difference is that the ion (electron) flow may strongly reduce the E x B drift-rate, 
to a level comparable with observationally deduced values. The time dependent shielding 
can result in a time variability of the observed drift-rate, including the apparent reversals 
of the drift direction. The latter effect can occur if the natural sampling rate (once per 
pulsar period) is too slow with respect to the drifting subpulses variability and results in 
an aliasing phenomenon. In fact, the apparent reversals of subpulsc drift direction can be 
explained by small variations (a few percent of the mean value) of the drift rate, which 
cause the P3 value to fluctuate around the relevant Nyquist boundary. 

2. Critical temperatures 

The value of the charge density above the polar cap heated by discharge bombardment is 
limited by the co-rotational GJ value. Since the number density of iron ions or electrons 
in the neutron star crust is many orders of magnitude larger than corotational values 
above the surface, then a thermionic emission from the polar cap surface is not simply 
described by the usual condition e c w kT s , where e c is the cohesive energy and/or work 
function, T s is the actual surface temperature and k is the Boltzman constant. Below 
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we consider pulsars with positively (ions) and negatively (electrons) charged polar caps 
separately. 



2.1. Iron critical temperature T{ 

In neutron stars with positively charged ($7 • B < 0) polar caps the outflow of iron 
ions is limited by thermionic emission and determined by the surface-binding (cohesive) 
energy e c . Let us consider, following the results of CR80, a general case of a pulsar inner 
accelerator in the form of a charge depletion region rather than a pure vacuum gap. 
According to their Eq. (8) the outflow of iron ions can be described in the form 

-exp(3Q-£r), (1) 



PGJ V kT > 

where pi < po.j is the charge density of outflowing ions. At the critical temperature 

T - = ss- < 2 » 

the ion outflow reaches the maximum value p = pej (Eq. (A.l)) permitted by the force- 
free magnetospheric condition. The numerical coefficient equal to 30 in Eqs. 0} and (J2J) 
is determined from the tail of the exponential function with an accuracy of about 10%. 
Thus, for a given value of the cohesive energy e c , the critical temperature Tj is also 
estimated within an accuracy of about 10%. Different values of the cohesive energy e c 
obtained by different authors lead to different values of critical temperatures. According 
to calculations of Abrahams & Shapiro (1991) 

T, = (6x 10 5 )6°- 73 (PiP_i 5 ) - 36 K, (3) 
while calculations of Jones (1986) give values five times lower 

T t = (1.2 x 10 5 )6 - 7 (PiP_i5) - 36 K, (4) 

where the parameter b = B s / Bd is described in Appendix A.l. (see also Eqs. and 
in Gil & Melikidze 2002; GM02 hereafter). 

Below the critical temperature Ti the charge-depleted acceleration region will form, 
with the accelerating potential drop AV — r\AV max ^ where AV max is the pure vacuum 
gap potential drop (Eq. (A. 2)), and the shielding factor can be defined/expressed in the 
form 

r, = 1 - Pl /p GJ = 1 - exp[30{l - TJT S )}. (5) 



2.2. Electron critical temperature T e 

Let us now consider pulsars with negatively charged polar caps (ft ■ B > 0), called 
"antipulsars" by RS. It is conventionally assumed that in such case a stationary free flow 
of electrons with the corotational GJ density (Eq. (A.l)) exists. This is the so-called space 
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charge limited flow (SCLF), in which the accelerating potential drop arises due to the 
dipolar field line curvature and/or inertia (Arons & Sharleman ll979l Arons H^81(l . Such a 
free flow requires that the electron work function w is completely negligible. However, it 
is known that w is of the order of 1 keV (see Eq. (C.6)), which is comparable with the ion 
cohesive energy e c (Eq. (2)). It is therefore possible that the electron flow is determined 
mainly by thermoemission. If, similarly to the ions case, the electron flow is slightly 
below the GJ value, then the effective potential drop just above the polar cap will be 
dominated by a small depletion of negative charge. This can happen if the actual surface 
temperature T s is slightly smaller than the critical electron temperature T e (Eq. ©)■ In 
fact, as demonstrated by Usov & Melrose (1996; see also Appendix C in this paper), the 
flow of electrons ejected from the polar cap surface by thermionic emission provides the 
GJ charge density if T s > T e ~ 0.04w/fc. One can therefore write electron analogues of 
Eqs. J3J and © in the form 

T - = 5BP (6) 

and 

r) = 1 - Pe/PGJ = 1 ~ exp[25(l ~ T e /T s )}, (7) 

where p e is the charge density of thermionic electrons, T e is the critical surface temper- 
ature and T s < T e is the actual surface temperature. 

The critical electron temperature T e is determined in terms of basic pulsar parameters 
by Eq. (|G.8(1 in Appendix C. Using Eqs. (| A. 2|> and <|A.3() we can rewrite this equation in 
the form 

T e ~ (5.94 x lO 5 )^^ 016 ^ K. (8) 

As one can see, the values of the critical electron temperature T e are close to the 
critical ion temperature T% obtained by Abrahams & Shapiro (1991). The observational 
estimates of polar cap surface temperatures T s based on first results from XMM satelite 
indicate values above 10 6 K (Becker & Aschenbach 2002). Therefore, the parameter b in 
Eqs. ©, Q and (JSJ) must be considerably larger then unity, implying a strong non-dipolar 
surface magnetic field (Eqs. (A. 2) and (A.3)). 

The original RS vacuum gap model is known to have a fundamental problem with 
the cohesive (binding) energy e c , which is apparently too low to bind the ^Fe ions in the 
uppermost layer of the polar cap surface (for review see Usov & Melrose 1995). However, 
Gil & Mitra (2001; GM01 hereafter) and GM02 argued recently that RS-type vacuum 
gap can form above a positively charged polar cap, provided that the surface magnetic 
field B s = bBd (Eqs. (A. 2) and (A.3)) is very strong (about 10 13 G) and non-dipolar in 
nature (radius of curvature 1Z < 10 6 cm). In the partially shielded acceleration region 
(r/ < 1), the actual surface magnetic field can be a few times lower, although still much 
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higher than the conventional dipolar field Bd (Table 1). In fact, one can show that the 
minimum surface magnetic field necessary for the formation of a pure vacuum gap (with 
r] = 1) obtained by GM02 (their Eqs. (8) and (14) and Fig. 1) should now be multiplied 
by a factor r? 57 . For six pulsars listed in Table 1, the average value of this factor is about 
0.3 and the average value of the parameter b = B s /Bd is about 4 and the average value 
of B d = 3.3 x 10 12 G, while the average value of B s = 8.7 x 10 12 G. 



3. Thermostatic regulation of the actual surface temperature T s 

Let us consider a quasi-equilibrium state when cooling by radiation balances heating due 
to electron bombardment of the polar cap surface 

oT s 4 = r]j acc m e c 3 n GJ , (9) 

where the shielding factor r\ is determined by Eqs. (J5J or (7), the Lorentz factor 7 acc is 
determined by Eq. (A.12) 1 , iiqj = PGj/e = 1-4 x 10 11 6(F_i 5 /P) - 5 cm" 3 , m e is the elec- 
tron mass, e is the elementary charge, c is the speed of light and er is the Stefan-Boltzman 
constant. It is straightforward to obtain the expression for the quasi-equilibrium surface 
temperature in the form 

T s = 4.34 x 10 6 P~ 1/7 Kl /7 ri 3/7 K. (10) 

Inverting this equation we obtain the expression for the shielding factor rj in terms of the 
surface temperature T s 

r) = 3.25 x 10- 2 Tq /3 pI /3 Kq 2/3 , (11) 

where T@ = T s /10 6 K. This expression describes the shielding factor rj, in terms of the 
balance Eq. © independently of the general definition (Eqs. or 10). Since Eq. (|11|) 
and Eq. J5J) or Eq. have to be satisfied simultaneously, then the actual surface tem- 
perature T s will be thermostatically self-regulated within a narrow range around the 
quasi-equilibrium value (Eq. In fact, a slight decrease of T s in Eq. © causes an 

increase of the shielding factor r/ (due to a smaller number of thermionic ions). This in 
turn causes an increase of T s in Eq. (16), due to a less shielded accelerating potential 
drop and more intense heating by discharge bombardment. 

Figure 1 shows values of the shielding factor r\ calculated from Eqs. (JSJ represented 
by two dashed lines, and l[ll[) . represented by six numbered solid lines. We introduced 
a parameter S = — T§ shown on the horizontal axis, which describes a difference 
between the iron critical temperature T{ = x 10 6 K and the actual surface temperature 



1 Although the Lorentz factor 7 acc was calculated within a framework of the NTVG-ICS model 
of GM01 (see also GM02 and Appendix A), its value is not strongly dependent on the adopted 
model of acceleration region (therefore T s is not strongly dependent on the particular inner gap 
acceleration model). 
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Fig. 1. Potential shielding factor -q calculated from Eq. JSJ - two dashed lines, and 
from Eq. (|I1|) - six solid lines numbered from 1 to 6, versus the temperature difference 
6 = (Ti — T s )/10 6 K, where 7j is the ion critical temperature (Eq. ©) and T s < % is 
the actual surface temperature. 

T s = Tq x 10 6 K. Using this parameter we can rewrite Eq. JHJ) in the form r\ = 1 — 
exp [— 308/ (TiQ — 6)), and two dashed lines presented in Fig. 1 correspond to limiting 
values TiQ =1 and 2, respectively 2 . Similarly, Eq. Ijlljl can be rewritten in the form 
77 = 3.25 x 10 _2 (7i6 — 5) 7 ^TZ e 2 ^ P 1 1 ^ 3 , and six numbered lines presented in Fig. 1 
correspond to different (indicated) combinations of and TZq values (Pi = 1 was used 
in all cases). The actual values of ?y correspond to the shadowed region limited by the 
two dashed lines and the two numbered lines 1 and 6 from the top and the bottom, 
respectively. If the iron critical temperature T, is not much lower than 10 6 K and ©) 
and TZq is not much smaller than 0.1, then the value of 77 should lie in the range 0.03-0.8. 
This suggests that the actual value of the parameter S should range between a fraction of 
a percent to several percent at most. We therefore conclude that a difference between Ti 
and T s (or T e and T s ) should be of the order of 10 4 K. In the next section we constrain 
both Ti (or T e ) and T s from the observationally deduced drift-rates for a number of 
pulsars. The inferred surface temperature T s and the difference between Ti (or T e ) and 
T s agree very well with a general estimate obtained here. 

2 Our conclusions are not very sensitive to these limiting values, provided that they are not 
much different from 1 and 2 (see Sect. 4). 
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The model of a charge-depleted acceleration region described by Eqs. (JIJ and Eq. (J3J) 
for ions and Eq. (|7J) for electrons is highly sensitive to both T^ e and T s temperatures, 
where T he means either T t (Eqs. ©, QJ) or T e (Eq. ©). The inequality T s < T ije is often 
used as a thermal condition for the "vacuum gap formation" (GM00; GM02; Abrahams & 
ShapiroEHH Usov & Melrose lT§551ITMl)|l and its exact meaning is worth understanding. 
When the surface temperature T s is only 10% below Xi. e , then 77 ss 0.95, that is pi, e is 
only a few percent of the corotational GJ charge density and the situation is close to the 
pure vacuum case of RS (see discussion below Eq. (8) in CR80 for details). If this is the 
case, then the potential drop AV^ « AV ma x is being screened mostly due to the intense 
production of e _ e + pairs, that is, the total charge density p t < pqj (equality holds when 
a spark plasma filament is fully developed), where p t = pi + p+ or p t = p e + p_, and p + 
and p- are charge densities of positrons and electrons produced in e~e + pairs in spark 
discharges, respectively. The back flow of high density electrons or positrons accelerated 
to relativistic energies by the vacuum potential drop, heats strongly the surface beneath 
the sparks. As a result, the surface temperature T s increases, which may intensify the 
thermionic ion/electron outflow, additionally screening the potential drop to the level 
AV = rj AVmax- At this stage, the cooling process should slightly dominate, because 
of weaker heating due to the back flow of lower density p± = rjpcj and lower energy 
7 = 7]AV max /m e c 2 electrons/positrons. Thus, the actual surface temperature should be 
thermostatically self-regulated. If heating and cooling are nearly in balance, then the 
quasi-equilibrium surface temperature T s is very close to Ti (Fig. 1), with the difference 
being less than about 1%. Thus, contrary to the suggestion of CR80 (see discussion below 
their Eq. (8)), the conditions close to a pure vacuum gap will never be established in the 
quasi-equilibrium stage. Nevertheless, the potential drop arising due to a small charge 
depletion above the polar cap is high enough for the non-stationary sparking scenario to 
be realized in the way similar to that envisioned by RS. 

The sparking discharge terminates if the total charge density pt = Pi, e + P± reaches 
the GJ value (Eq. (1)), where pt. e is described by Eqs. (5) and 10, respectively, and 
time dependent p±(t) is the charge density of produced electrons/positrons. A slight 
adjustment to pi t£ is also expected in the final stage of a spark development due to 
the surface temperature variations. More exactly, the discharge terminates when the 
screened potential drop AV = (1 — pt/ PGj)AV max = (t) — p±/ PG,j)AV max falls below 
the threshold value AV mm = j m in(mc 2 /e) « 2.5 x 10 3 T 6 1/3 ^ 1 / 3 (mc 2 /e) w 5 x 10 8 V, 
inhibiting further e~e + pair production (see Appendix A.l. for determination of AV max 
and j min ). Since AV max = AV ICS - (10 11 - 10 12 ) V (see Eq. (11) in GM02 or Eq. (A.6) 
in this paper), then AV m i n / AV max ~ 10 -3 ~ (77 — P±/pgj)- Therefore, the discharge 
terminates if p± = r/pcj, which can be a small fraction of GJ density. In fact, the 
condition 7 acc > 7 m i„ leads to a lower limit r\ > 10~ 4 bTZg 1 TQ : where 7 acc is determined 
in Eq. (A. 15). Since the product blZ^ 1 , describing small scale anomalies of the surface 
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magnetic field should be about 10, the absolute lower limit for 77 is about 0.001. It should 
be obvious from the above considerations that the sparking discharge can proceed in a 
similar way to that envisioned by RS, even if the parameter 77 is very low, say below 0.01. 
In reality, however, one should expect that typically 77 is of the order of 0.1 (Fig. 1) and 
thus pi te ~ 0.9pcj- Thus, one can say that the sparking discharge begins at the partially 
shielded potential drop AV ~ 0.1 AV m ax ~ 10 10 — 10 11 V and terminates when this 
potential drop is largely reduced by about two orders of magnitude (due to the screening 
by cascading production of e~e + pairs). 

At the final stage of a spark development, when heating is most intense, the maximum 
polar cap temperature that can be achieved is T s = T^ e and the spark discharge is ter- 
minated. Then, in the absence of heating from a pair production mechanism, the surface 
temperature should drop fast enough to enable the next discharge to recur within a suf- 
ficiently short amount of time. The e-folding cooling time r coo / is estimated in Appendix 
B.2. as being of the order of 1 /xs (Eq. B.14). Thus, the surface temperature can drop by 
a factor of e in a time interval of the order of ~ ps. However, in our scenario it is only re- 
quired that the temperature drops by a few to several percent, and thus a corresponding 
cooling time is at least 10 times shorter. This is comparable with a gap emptying time (or 
transit time) t ~ h/c, which is of the order of 100 ns or less (e.g. Asseo & Mehkidze ll998l 
GS00). On the other hand, the heating time scale due to electron/positron bombardment 
is about (10 — 40)h/c (RS, see also Appendix B.3.). Therefore, sparking discharges should 
be easily able to cycle on and off on a natural thermal timescales, comparable to those 
described by RS model of a spark development (exponentiation of charge density from 
p « to p = p GJ ). 



4. Subpulse drift in pulsars 

We now compare the observationally deduced drift-rates with theoretical predictions of 
the E x B drift model determined by Eqs. (A.13-A.15). The input parameters and the 
observationally deduced parameters are summarized in Table 1. The most important are: 
the drift periodicity P3 and the circulational (tertiary) periodicity P3, which are related 
to each other by the number of circulating sparks N = P3/P3. If necessary/possible 
these values are corrected for aliasing effects. In the first three cases listed in Table 1 we 
rely on a robust analysis of real observational data, while in the next three cases (PSRs 
2303+30, 2319+60 and 0031-07) we deduce drift parameters by reproducing the actual 
pulse sequences using the simulation model and fitting the free model parameters to the 
observed values (see GS00 for details). 
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4.1. PSR 0943+10 



This pulsar was extensively analyzed and interpreted recently by DR99, DR01, Asgekar & 
Deshpande (2001) and Gil & Sendyk (2003). These authors revealed clearly that drifting 
subpulses in this pulsar result from a system of N = 20 subbeams (sparks) circulating 
around the magnetic axis at the periphery of the pulsar beam (polar cap) . Each subbeam 
(spark) completes one full circulation in P3 = 37.35Pi, therefore P3 = P3/N = 1.87P%, 
where Pi is the pulsar period. These precise estimates 3 were possible due to the successful 
resolving of aliasing near the Nyquist frequency f% ~ 0.5/Pi (or P3 ~ 2Pi). 

At the periphery of the polar cap the circulation distance d ~ r p ~ h, where r p is 
the polar cap radius (Eq. (|A.4Jl ) and h is the height of the acceleration region (e.g. 
Eqs. I|A.8|) and (|A.9jl ). Thus equation l|A.13|l gives in this case P 3 = rj~ 1 P 1 {(r p /h) - 1]. 
In order to compare this equation with the observed value P3 = 37.35Pi, we have to 
estimate the complexity parameter r p /h, representing the ratio of the polar cap radius 
to the characteristic spark dimension V ~ h at the polar cap periphery (GS00). First, let 
as assume after RS that the typical distance between adjacent sparks at the periphery 
of the polar cap is also about h. If this is the case, then 2nd = 2n(r p — h) = 2Nh, 
which leads to r p /h — 1 = N/ir = 20/ir or r p /h rj 7.4 (see also GMM02b who obtained 
r p /h = 7.34 ± 0.8 for this pulsar; in consistence with approximate Eq. (11) in GS00, 
which gives r p /h — 6.8). Therefore, with r p /h = 7.4 we obtain P3 = QAPx/r) = 37.35Pi, 
which gives the value of the shielding parameter 77 = 0.17, that can be compared with 
theoretical predictions using Eq. (J5J or J7J for partial electron flow of ions or electrons, 
respectively. 

The actual quasi-equilibrium surface temperature is described by Eq. which for 
r] = 0.17 and TZ 6 = 1 gives T s = 1.646x 10 6 K. On the other hand, from Eqs. © and we 
obtain P/T, = 1.00625 and T e /T s = 1.0075, which gives critical iron T t = 1.656 x 10 6 K 
and electron T e = 1.658 K temperatures, respectively. Note that AT = T^— T s ~ T e —T s ~ 
10 4 K, in agreement with the general prediction made in Sect. 3. The comparison of Ti 
with Eq. J21 and Q for iron critical temperatures and T e with Eq. (JHJ for electron critical 
temperature implies that b = 2.8, 29.7 and 18 or B s ~ 10 13 G, B s = 11.2 x 10 13 G and 
B s = 7 x 10 13 G, respectively. As one can see, the calculations of the cohesive energy by 
Jones (1986 ) represented by Eq. |@} lead to values of B s largely exceeding the critical 
quantum field B q = 4.4 x 10 13 G, which is believed to represent the photon splitting 
limit above which the plasma necessary for generation of pulsar radio emission can not 

3 The precise values of P 3 = (37.35 ± 0.52)Pi and P 3 = (1.87 ± 0.026)Pi were given by 
DR99 and DR01 (at 430 and 111 MHz). Asgekar & Deshpande (2001) analyzing the 35-MHz 
observations of this pulsar confirmed that P3 ~ 37Pi. However, Rankin et al. (2003) analyzing 
103/40-MHz Pushchino observations demonstrated that both P3 and P3 values can slightly vary 
between different observing session, with differences being of the order of a few percent. 
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be produced (e.g. Zhang & Harding 2002 and references therein). On the other hand, 
calculations of Abrahams & Shapiro (1991) for the iron critical temperature represented 
by Eq. © give a quite suitable value of the surface magnetic field B s ~ 10 13 G (Table 
1). The electron case represented by Eq. © can be ruled out, since the corresponding 
B s = 7 x 10 13 G exceeds largely the photon splitting limit. 

4.2. PSR B0809+74 

This pulsar was recently analyzed by van Leeuwen et al. (20.QJ3, who argued that the 
number of circulating beams N > 15 and the circulation period P3 > 150Pi. Assuming 
again that the observed drifting subpulses in this pulsar correspond to the periphery of 
pulsar beam (polar cap) we have P3 = i]^ 1 Pi[(r / h) — 1] and r p /h — (N/ir) + 1> 5.8 (see 
Sect. 4.1). Thus, the shielding factor 77 > 0.032 and using Eq. J3J) we obtain Ti/T s = 1.001, 
while Eq. Q gives T e /T s = 1.0013. The actual surface temperature can be estimated 
from Eq. (TTJ, which for r\ = 0.032 and Tl 6 = 1 gives T s = 0.955 x 10 6 K. Therefore we 
have Ti = 0.956 x 10 6 K and T e = 0.957 x 10 6 K, respectively (notice that in this case 
AT = Ti - T s = T e - T s ~ 2 x 10 3 K). These values can now be compared with Eqs. 0, 
J3J and (JSJ), which give B s — 4x 10 12 G (less than the minimum value 0.1B q required by 
near-threshold conditions described in Appendix A.I.), B s = 3 x 10 13 G (too close to B q ) 
and B s = 7 x 10 12 G, respectively (OK). The latter value of the surface magnetic field is 
very suitable, suggesting that PSR 0809+74 is a good candidate for a drifting subpulse 
pulsar with a positively charged polar cap (antipulsar in the terminology of RS). 

4.3. PSR B0826-34 

The drifting subpulses in this remarkable pulsar were first observed by Biggs et al. ([1985JI . 
The average profile is very broad and consist of the main-pulse (MP) and the interpulse 
(IP). The single pulse emission occurs practically at every longitude of the 360 degrees 
pulse window, implying that this pulsar is an almost aligned rotator (i.e. the inclination 
angle a between the rotation and magnetic axes is very small). At 645 MHz Biggs et 
al. I|1985f) revealed 5 to 6 bands of drifting subpulses, swinging across the MP window. 
Individual subpulses show a variable drift rate, including changes of the apparent drift 
direction approximately every 100 pulses (see Fig. 2 in Biggs et al. 1985). Such patterns of 
drifting subpulses are inconsistent with the E x B drift model, unless the observed effects 
result from aliasing and small changes in the drift-rate causing P3 to fluctuate around 
the relevant Nyquist boundary (as in the case of PSR 2303+30; see Fig. 3 in GS00). 
Recently Gupta et al. (|2003J) attempted to resolve the possible aliasing in PSR 0826- 
34. They reobserved this pulsar at 318 MHz using the GMRT and obtained a sequence 
of 500 good quality single pulses. At this frequency 6 to 7 bands of drifting subpulses 
appeared, behaving similarly to those 5-6 bands observed by Biggs et al. (19831) a t a 
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higher frequency. The careful analysis revealed that subpulse-associated sparks circulate 
always in one direction consistent with the E x B drift model, but the observed patterns 
are determined by the aliasing corresponding to P 3 ~ Pi (thus around the fluctuation 
frequency / 3 ~ 1/Pi). 

The observed longitudinal drift rate Dq = dip/dt — P2/P3 corrected for aliasing (P3 ~ 
Pi) is Dq « P2/P1, where P2 « 26° is the typical distance between adjacent subpulses, 
where <p is the observed subpulse longitude. The circulational periodicity P3 = 360°/P/ r , 
where D r — dx/dt is the intrinsic azimuthal drift rate (see Appendix A. 2.). In the case 
of an almost aligned rotator (a ~ 0°) one can write D r ~ Dq = 26°/Pi and thus P3 w 
(360°/26°)Pi = 13.85Pi. Therefore, in PSR 0826-34 there are N = P 3 /P 3 w P3/P1 = 14 
sparks circulating around the pole at the rate D r rj 26° /Pi. 

The detailed analysis of single pulses performed by Gupta et al. H2003|) revealed also 
that sparks contributing to the MP emission of B0826-34 circulate at d ~ 0.5r p (the 
middle of the polar cap), while the IP emission occupies the periphery of the polar cap. 
Thus, according to Eq. I|A.13(1 with s = 0.5, the circulation period P3 = t]^ 1 Pi(0.5r p / 7i) 2 , 
which can be compared with the observationally deduced value P3 ~ 14Pi. To estimate 
the value of (r p /h) let us notice that in this case 2nd = 2NT>, where d ~ 0.5r p and 
D ~ 0.5/i is the characteristic spark dimension in this region of the polar cap 4 . Therefore 
r p /h = N/ir = 14/tt ~ 4.5, thus P3 = 5rj~ 1 Pi w 14Pi and, in consequence, the shielding 
factor r) = 0.357. Now using Eqs. © and we obtain T,/T s = 1.0147 and T e /T s = 
1.0188, respectively. The actual surface temperature can be estimated from Eq. i|10fl . 
which for TZq = 1 and 77 = 0.36 gives T s — 2.565 x 10 6 K. Consequently, we get Tj = 
2.469 x 10 6 K and T e = 2.475 x 10 6 K, respectively (note that AT = Tj — T s w T e - T s ~ 
3 x 10 4 K). Comparing these critical temperatures with Eqs. J3J, Q and JHJ we obtain 
P s = 1.4 x 10 13 G, 1.5 x 10 15 G and 5 x 10 13 G, respectively. Thus, again the calculations 
of the cohesive energy by Abrahams & Shapiro l)1991|l represented by Eq. © seem to be 
the most adequate, while the calculations of Jones (1986) lead to B s > 10 15 G. Also the 
electron case with B s = 5 x 10 13 G > B q can be excluded in this case. 

The observed time variability of the apparent drift rate, including the change of 
the apparent drift direction in PSR B0826-34 can be explained by time variations of 
the shielding factor 77 (Eqs. JSJ, Q, (A. 7) and (A. 8)), which, according to Eq. (fTTfl . 
implies a time variability of the surface temperature T s . The detailed analysis of sequences 
of drifting subpulses in PSR 0826-34 performed by Gupta et al. (2003) indicates the 
following variations over about 100 pulse sequences: the azimuth drift rate D r varies 
from 25°/Pi to 27°/P 1; drift velocity v d from 7.5 m/s to 9 m/s, P 3 from 1.03Pi to 
0.95Pi and P3 from 14.4Pi to 13.3Pi. The change of the apparent drift direction due to 
4 We assume here that T> ~ h at the polar cap boundary (d ~ r p ) but due to convergence of 
the planes of field lines (driving the spark avalanche) towards the pole T> ~ 0.5/i in the middle 
of the cap (d ~ 0.5p). 
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aliasing occurs when P3 = l.OPi and thus P3 = lAPi. These approximately 100 pulse 
cycles repeat in a quasi-periodic manner. During each cycle the shielding factor increases 
by a few percent (around rj = 0.36) and the surface temperature drops by about 2000 
K (around T s = 2.47 x 10 6 K). Notice that similar (though erratic) variations of these 
parameters, with an amplitude of about a few percent, were observed in PSR 0943+10 
(see footnote 3). 

4.4. PSR 2303+30 

GS00 reproduced the drifting subpulses in this pulsar (see their Fig. 3) and argued that 
in this case r p /h ~ 5, TV — 12 and P 3 w 23Pi, which is consistent with P 3 ~ 1.9Pi 
obtained by Sieber & Oster H1975f) . For the periphery of the polar cap we get 77 = 0.16, 
T s = 1.854 x 10 6 K, Ti = 1.864 x 10 6 K and T e = 1.865 x 10 6 K. Comparison with 
Eqs. ©, and © gives B s = 9.6 x 10 12 G (OK), B s = 8.7 x 10 13 G (too high) and 
B s = 3.7 x 10 13 G (probably too high), respectively. Thus again, the calculations of the 
cohesive energy by Abrahams & Shapiro (1991) seem to be most adequate. 

This pulsar also shows an apparent change of the subpulse drift direction (see Fig. 3 in 
GS00), which can be explained by small changes of: the azimuth drift-rate from 14.°4/Pi 
to 15.°5/Pi, the drift velocity Vd from ~ 10 m/s to ~ 11 m/s, P3 from 2.1Pi to 1.9Pi and 
P3 from 25Pi to 23Pi (notice, that similar variations of these parameters were observed 
in PSRs 0826-34 and 0943+10). The corresponding change of the shielding factor rj is 
about one percent, which implies the change of the surface temperature T s of about 
1000K (around T s = 1.85 x 10 6 K). 

4.5. PSR 2319+60 

GS00 reproduced the drifting subpulses in this pulsar (see their Fig. 4) and argued that 
in this case N = 9, P3 ~ 70Pi, P3 ~ 7.8Pi and r p /h ~ 5. For sparks drifting at the 
periphery of the polar cap we get 77 ~ 0.071, T s = 1.236 x 10 6 K, T, = 1.239 x 10 6 K 
and T e = 1.239 x 10 6 . Comparison with Eqs. ©, © and © gives B s - 5 x 10 12 G 
(this value is lower than Bd = 8 x 10 12 G and should be rejected), B s = 5 x 10 13 G 
(this value is higher than the photon splitting limit B q and should also be rejected) 
and B s — 1.36 x 10 13 G. Thus, only the latter case seams acceptable and therefore 
PSR 2319+60 should be considered as the pulsar with a negatively charged polar cap 
(antipulsar in RS terminology). 

4.6. PSR 0031-07 

GS00 reproduced drifting subpulses in this pulsar (see their Fig. 5) and argued that in 
this case iV > 5, P3 > 34Pi, P 3 ~ 6.8P1 and r p /h ~ 4. For the peripheral sparks we get 
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J) > 0.089, T s = 1.551 x 10 6 K, T l = 1.556 x 10 6 K and T e = 1.557 x 10 6 K. Comparison 
with Eqs. ©, Q and © gives B s = 7.2 x 10 12 G (OK), B s = 7.9 x 10 13 G (too high) and 
B s = 2.2 x 10 13 G (probably too high). Thus again, the calculations of cohesive energy of 
Abrahams & Shapiro (1991) seem to be most adequate and calculations by Jones (1986) 
are not suitable within our model. 

5. Discussion and conclusions 

Recent progress in the analysis of drifting subpulses in a number of pulsars (DR99, DR01, 
Asgekar & Deshpande 2001, van Leeuwen et al. 2002, GMMOlb, Gil & Melikidze 2002 
and Gil & Sendyk 2003) provided a strong support to the canonical RS pulsar model, 
which relates the spark plasma filaments with subpulse-associated subbeams of radio 
emission. The observed subpulse drift is naturally explained by the inevitable E x B 
drift of these plasma filaments, at least qualitatively. However, on the quantitative level, 
the E x B drift in a pure vacuum gap is too fast as compared with observations. In 
this paper, we consider a general concept of the charge depletion region (rather than a 
pure vacuum gap), with a thermal outflow of ^Fe ions, coexisting with the generation 
of e~e + plasmas (CR80). The presence of an ion flow decreases the accelerating gap 
potential drop and, in consequence, the E x B drift-rate, as well as the amount of surface 
heating due to back-flowing relativistic electrons (positrons). As originally demonstrated 
by CR80, the extreme sensitivity of the potential drop to the surface temperature makes 
the gap thermostatically self-regulating, whenever there is both an ion outflow and e~e + 
discharge plasma. We argued that when heating and cooling are nearly in balance, the 
surface temperature T s oscillates within a very narrow range around its quasi-equilibrium 
value, slightly below the critical temperature T{ (T e ) above which the ion outflow reaches 
the co-rotational GJ density. 

Our general model of the charge depleted acceleration region can be applied not 
only to pulsars with positively charged polar caps (SI ■ B < 0), but also to pulsars 
with negatively charged polar caps (SI ■ B > 0), which are usually interpreted within 
the so-called stationary flow models (e.g. Arons 1981 and references therein) in which 
electrons flow freely from the polar cap surface at the GJ charge density. Within the 
model of a partially shielded acceleration region the ^t^e ions are only marginally bound 
at a temperature T s <Ti, but, nevertheless, the non-stationary sparking scenario can be 
realized in a way similar to that envisioned by RS in their pure vacuum gap model 
(with strong binding assumed). We suggest that such marginally bound ions are not 
much different from marginally bound electrons (at T s ~ T e ). In fact, RS-type non- 
stationarity (sparks) can also arise if the thermionic electron flow from the surface is not 
exactly at the GJ corotational charge density. As we demonstrated, even small departures 
from GJ charge density at the surface temperature T s w T e can result in the sparking 
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discharge of the shielded accelerating potential drop above the polar cap. Thus we suggest 
that both pulsars with positively and negatively charged polar caps can develop non- 
stationary inner acceleration regions, similar to those invoked by RS. In consequence, 
the spark-associated coherent radio emission due to instabilities in the non-stationary 
e~e + secondary plasma (e.g. Usov 1987, Asseo & Melikidze 1991, Melikidze et al. 2000) 
may originate similarly in both pulsars (fJ • B < 0) and "antipulsars" (O • B > 0). 

We estimated values of the polar cap surface temperature for six pulsars (Table 1). 
The estimated values of T s lie in the range between ~ 1.0 to ~ 2.5 million K. It is worth 
noting that in a pure vacuum gap (in RS with r\ = 1), the values of T s would be higher 
by a factor of about 2, thus in some pulsars from our list exceeding 4 x 10 6 K at the polar 
cap. Present and future x-ray satellite spectral measurements should be able to constrain 
maximum surface temperatures of the polar cap heated by intense spark discharges. One 
should also mention that our estimates of T s were obtained under canonical assumptions 
that TZq — 1 (e.g. RS). However, if the radius of curvature of surface magnetic field lines 
is much smaller than 10 6 cm, then the values of T s can be lowe even by a factor 2 — 3. 
Thus, our estimates of T s presented in Table 1 represent the upper limits. 

Our method is capable of determing both the actual surface temperature T s and crit- 
ical temperature Ti or T e , above which the thermionic emission of iron ions or electrons 
reaches the corotational GJ value. We found that radio pulsars operate at T s approxi- 
mately 10 K lower than Ti or T e . Comparing our observationally deduced values of crit- 
ical temperatures with estimates based on the calculations of the cohesive energy and/or 
work function, we found the required values of the surface magnetic field B s = bBd (Table 
1). Within our model we can exclude the calculations of the cohesive energy by Jones 
(1986 ), represented in our paper by Eq. Q. In most cases calculations of Abrahams & 
Shapiro (1991) represented by Eq. Q give suitable results, except PSR 0031-07 which is 
the best candidate for drifting subpulse pulsars with a negatively charged polar cap. Also 
PSR 0869+74 seems to belong to this category. Interestingly, these two pulsars require 
the lowest shielding factors rj (Table 1) to fit the observed and the theoretical drift-rates. 

The sparks rotate due to the E x B drift in the same direction as the neutron star 
rotates (lagging behind the polar cap surface). As a matter of fact, the particle drift in the 
inertial frame cannot exceed the stellar rotation speed within any reasonable model. Thus, 
the observed direction of the subpulse drift is fixed for a given pulsar and depends only 
on whether the line-of-sight trajectory passes poleward or equatorward of the magnetic 
axis. However, two pulsars from our sample demonstrate an apparent change of the drift 
direction. We argued that sense-reversals can be explained by small changes of the drift- 
rate causing the P3 value to fluctuate around the relevant Nyquist boundary: P3 rs 1P\ 
in the case of PSR 0826-34 and P 3 ~ 2P 1 in case of PSR 2303+30, corresponding to 
fluctuation frequencies equal to 1 cycle/Pi and 0.5 cycle/Pi, respectively. We believe 
that these changes are due to small variations of the polar cap surface temperature. 
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PSR 0540+23 probably also belongs to the category of pulsars with P3 w 1P 1; showing 
apparent reversals of the sense of subpulse drift (Nowakowski 1991), although its drifting 
patterns are much more chaotic than in PSR 0826-34. 

Up to now, it was believed that in conal profiles (Rankin 1983) the values of P3 
periods covered the range from about 2 to about 15 pulsar periods Pi (corresponding to 
fluctuation frequencies between about 0.5/Pi (Nyquist frequency) to about 0.07/Pi). It 
is well illustrated in Fig. 4 and Table 2 in Rankin (1986). In this paper we broadened this 
range by adding at least one pulsar PSR 0826+23 with P 3 ~ P 1 (Table 1), and probably 
PSR 0540+23 (Nowakowski 1991) also belongs to this category. Perhaps there are many 
more pulsars with a fast drift corresponding to P3 values well below "Nyquist limit" . 

It is worth emphasizing again that the subpulse drift is widely considered as a cru- 
cial clue towards solving the long standing mystery of pulsar radio emission. Despite 
its potential importance, this phenomenon has not attracted much theoretical attention 
beyond that of RS model. A simple explanation of this fact is that probably it is very 
difficult to propose a theory that would be as natural as the RS model, at least quali- 
tatively. Moreover, this is the only model offering quantitative predictions that can be 
compared with observations. We review in this paper a number of problems appearing 
when the existing theory is confronted with the current pulsar data and argue that all of 
them disappear when the original RS model is modified to include a thermionic outflow 
of ions or electrons, coexisting with the electron-positron plasma produced in spark dis- 
charges. We also emphasize the importance of the strong non-dipolar surface magnetic 
field driving the spark avalanches, as well as resonant inverse Compton scattering as the 
mechanism providing seed photons for these avalanches. 

Finally, we believe that the phenomenon of drifting subpulses is a manifestation of 
a more general phenomena related to sparking discharges of the ultra-high potential 
drop above the polar cap of radio pulsars. Therefore, the result of this paper strongly 
suggest that spark-associated plasma instabilities (Usov 1987, Asseo & Melikidze 1998) 
play an important role in generation of the observed coherent pulsar radio emission (e.g. 
Melikidze et al. 2000). 



Acknowledgements. This paper is supported in part by the KBN grant 2 P03D 008 19 of the 
Polish State Committee for Scientific Research. JAG acknowledges the renewal of the Alexander 
von Humboldt Fellowship and thanks G. Riidiger for hospitality at the Astrophysikalisches 
Institut Potsdam, where this work was done. We also thank E. Gil and U. Maciejewska for 
technical help. 



Gil, Melikidze, & Geppert: Drifting subpulses... 19 
Appendix A: Inner acceleration region 

A.l. Charge depleted acceleration region above the polar cap 

The acceleration potential drop above the polar cap results from the deviation of a local 
charge density p from the corotational GJ density 

= ^ 

where f2 is the pulsar spin vector and B(r) is the pulsar magnetic field, which should 
be purely dipolar at altitudes of a few stellar radii P* = 10 6 cm, where the pulsar radio 
emission is expected to originate (e.g. Kijak & Gil 1998 and references therein), but could 
be highly nondipolar (i.e. having a pronounced small-scale spatial structure) at the polar 
cap surface. For reasons of generality we describe the surface magnetic field by 

B s = bB dl (A.2) 

where 

B d = 2 x 10 12 (P!P_ 15 ) - 5 G (A.3) 

is the dipolar component of the pulsar magnetic field at the pole, P is the pulsar period 
in seconds and P_i5 = P/10~ 15 and b > 1. This defines the actual radius of the polar 
cap 

r p = 6-°- 5 10 4 Pf 5 cm, (A.4) 

(e.g. GSOO). We also introduce the curvature radius of nondipolar surface magnetic field 
lines TZe = P/P*< 1, where P* = 10 6 cm is the neutron star radius (note that for 
a purely dipolar field Rq ~ 300 near the polar cap). Within our model of the charge 
depleted region the local charge density p = (1 — rfpoj, where r) = 1 — pt/pcj ( or 
77 = 1 — Pe/ Pgj) and Pi (or p e ) is the charge density of outflowing thermionic ions (or 
electrons) . Within the acceleration region the potential V and the electric field E (parallel 
to the magnetic field B s ) are determined by one-dimensional Poisson equation d 2 V/dz 2 = 
—4tt(p — pgj) = ^VPgj, with the boundary condition dV/dz — at z = h = h max and 
dV/dz = E max at z = (see below for determination of h). Since pa.j ~ ±P s /cP, where 
the sign +(— ) corresponds to the ion (electron) flow, the solution of Poisson equation 
gives E{z) = ±T ] (4Tr/cP 1 )B s (h - z) and AV = f£ E(z)dz = ±t]{2ti / 'cP^B.h 2 or 

AV = ± V AV max , (A.5) 

where 

2tt 

AV max = —B s h 2 (A.6) 
cPi 

is the acceleration potential drop within a pure vacuum gap (RS) corresponding to 77 = 1 
{pi = or p e = 0). 



20 Gil, Melikidze, & Geppert: Drifting subpulses... 

The height h of the acceleration region with ultra-high potential drop AV described 
by Eqs. (A.l) and (A. 2) is determined by a quasi-steady discharge via magnetic conver- 
sion of high energy fko > 2mc 2 photons into e~e + pairs that develop cascading sparks. 
Thus h = l p h, where l p h is the mean free path of an energetic photon to produce a pair. At 
least several possible models of the inner acceleration region were considered by different 
authors (see Zhang et al. 2000; ZHM00 hereafter and GM02 for review). Among them, 
the most promising one with respect to the possibility of creating an effective accelera- 
tion region ("vacuum gap") seems to be the Near- Threshold (fko = 2mc 2 lZ/l p h) model, 
involving the resonant ICS seed photons (Tko = 2"fheB s /mc) in the strong magnetic field 
B s > 0.lB q ~ 5 x 10 12 G (see GM01 and GM02 for details of the so-called ICS-NTVG 
model) . Unlike in the RS case with curvature seed photons, the mean free path of ICS 
seed photons in ICS-NTVG model is approximately equal to the mean electron path 
l e = 0.00277 2 £? 12 1 T f r 1 to emit an energetic photon with energy fko — 2'yheB/mc ex- 
ceeding a pair formation threshold (ZHM00, GM01, GM02), where B 12 = B S /1Q 12 G. 
Since 2mc 2 lZ/h = 2^/fieB s / mc, then the minimum Lorentz factors 7 required for the ICS 
dominated pair production under the Near Threshold conditions is 

lmin = (2.18 x lO^TZeh-^PP^)- - 5 , (A.7) 

(see Eq. (2) and GM01 and references therein for details). Now using hies — h — l e with 
7 = 7mm we can obtain the expression for the height of the ICS dominated acceleration 
region with a strong magnetic field 

hjes = (8.6 x ltf)nl'h-\P ■ P^y^T- 113 cm (A.8) 

and thus 7min « 2.5 x 10 3 T 6 1/3 ^g /3 , which reproduces Eq. (10) in GM02 

hies = (5 x lO 3 )^- 5 ^-^- - 36 ^ cm, (A.9) 

if T 6 = T s /10 6 K is expressed by their Eq. (12) and 77 = £ = k = 1. According to Eqs. (2), 
(A.l) and (A. 2), the acceleration potential drop is 

AV=(8x lO 12 )^- 1 ^ 73 ^- 27 ^- 3 / 2 ^ 1 / 2 V, (A.10) 
which reproduces Eq. (11) in GM02 

AV ICS = (5.2 x io 12 )7zJ- u r 1 p- 1 - 22 pri°B 8 V, ( A - n ) 

again under conditions mentioned above. The electrons (positrons) can be accelerated by 
this potential drop to the Lorentz factor values 

_ £^ _ 

1.8 x lO^ft-^^'p- 3 / 2 ^ 1 / 2 ^- 273 . (A.12) 

It should be noted here that in the case of positively charged polar cap the full 
accelerating potential drop AV = rjAV max + (1 — f])AVscLF, where AVsclf ~ 3 x 



Gil, Melikidze, & Geppert: Drifting subpulses... 21 

10 11 P 1 _1 ^ 2 (i? s /10 12 G) 1 / 2 (ft-/10 4 cm) V is the potential drop arising due to ions inertia 
in the space-charge-limited flow (e.g. CR80). If h ~ hies < 10 3 cm (Eq. (A. 9)) and 
AV m ax ~ AV7cs (Eq. (A.ll)) then for 77 - 0.1 (Table 1) AV ~ 77AF TOax . Thus, within 
the acceleration region the SCLF potential drop can be neglected, even when the general 
relativistic (GR) effect of inertial frame dragging (Muslimov & Tsygan 1992) is taken into 
account. In fact, the GR- induced electric field Eqh, grows quasi-exponentially from zero 
at the surface to the maximum value at the height approximately equal to the polar cap 
radius r p 3> h < 10 3 cm. Thus, the value of the integral EgrcIz must be negliglible 
compared with AV (e.g. Eq. (A. 10)). 



A. 2. E x B drift in the acceleration region 

If the actual charge density p < pej (thus 77 < 1), then the plasma within the acceleration 
region does not corotate with the neutrons star. This is the so-called E x B drift, which 
results in a slow motion of e~e + plasma in the direction perpendicular to the planes of the 
local surface magnetic held lines. Effectively, any filamentary lateral structures (sparks) 
circulate slowly around the local magnetic pole with a velocity Vd — cAE±/ B s , where 
AE± is the component of the electric field caused by charge depletion Ap = po.j — p = 
rjpcj, B s is the local surface magnetic field and c is the speed of light. According to 
Eq. (30) in RS, this electric field AE± = AV/(r p — d), where AV is the potential drop 
described by Eqs. (A.l) and (A. 2), r p is the radius of the polar cap (Eq. (|A.4|l ). and 
d < r p is the circulation distance of sparks from the local magnetic pole. 

Each spark completes one full circulation around the pole in a time interval P3 ps 
2nd/vd seconds. Using Eqs. (A. 5) and (A. 6) we obtain the so-called tertiary periodicity 
(expressed in pulsars periods Pi) 

where s — d/r p . Consequently, the azimuthal drift-rate D r = dx/dt (where \ 1S the 
magnetic azimuth angle of the circulating spark) is 
360^3^ ^ 

We assume reasonably that the biggest contribution to the E x B drift occurs at the 
beginning of the spark discharge, when AV ~ r]AV max (Eqs. (A. 5) and (A. 6)) that is 
when the screening due to exponentially growing e~e + pair plasma is weak. This justifies 
the form of our Eqs. (A. 13) and (A. 14). 

Assuming d = r p /2 (thus s — 0.5), we have P3 = 0.25(r p /h) 2 , which for rj — 1 (a pure 
vacuum gap) reproduces Eq. (32) in RS to within a factor of 0.25. This discrepancy follows 
from the factor 2 obviously missing in their Eq. (31). In fact, RS used the approximation 
< ae >~ r p /2 expressed explicitely just above their Eq. (31), in which, however, they 
put d =< ae >= r p , inconsistent with their assumption. 
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As originally argued by RS, the spark model very strongly suggests that 

P 3 =NP 3 , (A. 15) 

where N is the number of circulating sparks, and P 3 = l// 3 is the basic drifting periodic- 
ity determined either visually or measured from the dominant feature / 3 in the fluctuation 
spectrum (e.g. Backer 1973). If necessary, the value of P 3 has to be corrected for aliasing 
to allow to describe adequately the relationship between the tertiary periodicity P 3 and 
the number of circulating sparks N. 

Appendix B: Heat transport at the polar cap surface 

B.l. Polar cap physics 

An inherent property of the polar cap is the presence of the ultra-strong magnetic field 
(ignored by CR80), which is almost perpendicular to the surface. The presence of such 
strong magnetic field affects the equation of state, as well as the transport processes in the 
outermost surface layers of the neutron star. The state of matter and the heat transport 
regime in the polar cap region is determined by a number of quantities described in the 
following subsections B.l.l. - B.l. 8: 

B.l.l. Surface density 

The existence of a magnetic field in the order of magnitude of ~ 10 13 G enhances the 
binding energy of the electrons. This changes drastically the density profile in the very 
surface layers, the pressure (in the Thomas-Fermi approximation) vanishes at the so- 
called zero pressure condensation density p s which may be regarded as the surface density 
for a 3-dim bec crystal (Rognvaldsson et al. 1993). The most recent result for that density 
is obtained by Lai (2001) 

Ps ~ 561AZ~ 3 ^B% 5 g cm" 3 . (B.l) 

Here we consider the neutron star surface to be composed of iron, i.e. A = 56 and 
Z = 26. Thus, the surface matter density is p s (-Bi2 = 1) = 4.45 x 10 3 g cm~ 3 and 
Ps(Bi2 = 10) = 7.05 x 10 4 g cm" 3 , respectively. 

B.1.2. State of aggregation 

The state of aggregation of the polar cap matter is determined by the Coulomb parameter 
r = (Ze) 2 /(akT), which measures the ratio of electrostatic and thermal energy; a is the 
Wigncr-Scitz cell radius oc p -1 / 3 . Once the Coulomb parameter exceeds T = 170, the 
ions form a crystal (see e.g. Slattery et al. 1980), the corresponding melting density is 



170 kT\ 3 3Am u 



{Ze) 2 ) 4tt 



= 76 Tg gem' 3 , (B.2) 
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where m u is the atomic mass unit. Since p m <C p s {B\2 = 1) up to temperature Tq 
the polar cap matter will not be melted but consist of crystallized iron. 
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= 3, 



B.1.3. Degree of ionization 

For p > 10 5 g cm' 3 the atoms are fully ionized due to the pressure independently of 
the temperature, because the mean volume of a free electron becomes smaller than the 
volume of the Is— orbitals. For p s (B\2 = 1 or 10) the mean ionization per iron atom in 
the temperature range Tq > 1 is 24 or 25, respectively (Schaaf 1988), i.e. the atoms in 
the polar cap are almost completely ionized. The iron ions are non-degenerated and not 
affected by the magnetic field, i.e. the phonon spectrum of the crystal remains the same 
as for the zero magnetic field (B = 0) case. 

B.1.4. Degree of degeneracy 

The Fermi temperature for the iron polar cap is given by 



(see e.g. Hernquist 1984). In the presence of a quantizing magnetic field the electron 
chemical potential is reduced. If only the lowest Landau level is populated, the Fermi 
temperature is 



Thus, for the typical densities and magnetic field strengths in the polar cap region the 
electrons there are in the state of a complete degeneracy, provided that the surface is not 
hotter than a few times 10 6 K. 

B.1.5. Degree of being relativistic 

That state of the electron gas is determined by the relativistic parameter x — pf(B = 
0)/m e c (Yakovlev & Kaminker 1994). The Fermi momentum at B = is given by 



For p s {B\2 = 10) we find pp ~ 8.8 x 10 18 g cm s 1 which is much less than m c. 
Therefore, the electrons in the polar cap region are non-relativistic. 

B.1.6. Number of populated Landau levels 

As shown by Yakovlev & Kaminker (1994) for conditions realized at the polar cap, the 
strong magnetic field > 10 12 G has a quantizing effect on the electron motion. How many 




(B.3) 




(B.4) 




(B.5) 
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Landau levels arc typically populated? The maximum number of these levels is given by 
(see e.g. Fushiki et al. 1989, Hernquist 1984) 

p*(B = 0)c p* /3 

n ™* ' 2ehB ~ U ^7- (B ' 6) 

For B12 « 10 and p < p s (Bi2 = 10) = 7 x 10 4 g cm -3 , only the lowest Landau level 
(n = 0) is populated. 

B.1.7. Debye temperature 

For temperatures larger than the Debye temperature 6d the specific heat of the bec 
crystal lattice can be approximated by its classical value Ci at tice = 3fcn;, where m is the 
number density of ions. According to Yakovlev & Urpin (1980), the Debye temperature 
as a function of the ion plasma frequency lu p is given by 

6 D = 0.45^ w 1.6 x 10 6 pl /2 K. (B.7) 

The presence of a strong magnetic field does not change 0d drastically. Since p s {B\2 = 
10) <C p%, we find for the polar cap layer T > Od and the specific heat of the crystal 
there can be regarded as caused by classical lattice vibrations. 



B.1.8. Depth of heat deposition 

The depth of the heat deposition due to the bombardment of the polar cap by ultra- 
relativistic electrons/positrons is measured in so-called radiation lengths (see CR80). 
For IgFe ions the radiation length y = 14 g cm~ 2 . Therefore the corresponding depth 
L w y/p s , and for B\ 2 — 10 we have L~2x 10~ 4 cm. 

B.2. Cooling timescales 

The polar cap surface is heated by the bombardment of back-flowing relativistic electrons 
produced in the pair producing spark discharge of the accelerating potential drop. Once 
the total charge density pt — Pi, e + P± approaches the GJ value, the intense heating 
ceases until the next spark recurs. In the absence of heating the region beneath the 
spark cools down rapidly. The cooling process was considered by CR80, who found the 
characteristic cooling time r c = 0.32(T S /10 6 K )~ 6 (/?™/10 5 gcm~ 3 ) 7 / 3 seconds, where T s 
is the actual surface temperature and p™ is the actual matter density just beneath the 
polar cap surface. Adopting T s = 10 6 K and p™ = 10 5 g cm~ 3 CR80 obtained a small 
fraction of a second for r c (which could be a very tiny fraction of a second if T s is slightly 
higher and p™ is lower). It seems therefore appropriate to reconsider their derivation, 
since much more precise estimates for the transport coefficients and their magnetic field 
dependence in the uppermost surface layer are currently available. In what follows we 
omit the superscript m denoting the matter density p m . 
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As we have shown in Appendix B.l. the magnetized (B\2 = 10) surface of the polar 
cap in isolated (non-accreting) neutron stars has a density p s rts 7 x 10 4 g cm -3 and at 
temperatures of T s > 10 6 K consists of almost completely ionized iron ions which form 
a solid bec crystal lattice having a Debye temperture well below 10 6 K. The electrons in 
that lattice form a degenerated non-relativistic gas which populates the ground Landau 
level only. 

Let us now estimate the characteristic cooling time for the strongly magnetized out- 
ermost surface layer of the polar cap. The cooling time can be deduced from the heat 
transport equation. The heat transport in the polar cap region can be considered in a 
very good approximation as purely parallel to the magnetic field lines, which are almost 
perpendicular to the surface. Therefore, the heat transport is well described by a one- 
dimensional heat transport equation with the use of the time coordinate t and the spatial 
length coordinate I (parallel to the magnetic field lines). Without sinks and sources of 
energy and together with the boundary condition that the heat is radiated away from 

9T 
dl 

the heat transport equation reads 



the very surface according to — k 4j£ = a T 4 (where a is the Stefan-Boltzman constant), 



BT d ( dT\ 

c m=di{ K m)- (R8) 

By assuming in the very thin surface layer an almost uniform heat conductivity, and 
using dT/dt « T/r c and d 2 T/dl 2 « T/L 2 , we obtain the e-folding cooling time 

Tcool = (^p) , (B.9) 

where L is the depth of heat deposition, k is the thermal conductivity and C is the 
specific heat per unit volume. 

Heat can be stored both in lattice vibrations and in the electron gas. In the parameter 
range under consideration these contributions to the specific heat are comparable, thus 

C = Clattico + Col • (B.10) 

For the number of electrons per nuclcon corresponding to ^Fe ions, the contribution of 
the non-relativistic degenerated electron gas to the specific heat is 

C e i « 1.05 x 10 U P6 /3 T 6 erg K _1 cm _3 , (B.ll) 

(see Landau & Lifshitz 1969). Since at the polar cap T > 9d the lattice specific heat is 
given by 

3k 

Catticc = 3 Hi k = — p w 

Am u 

AA x 10 12 /?6 erg K _1 cm- 3 . (B.12) 
Thus, we obtain for the specific heat per unit volume at the polar cap region 
C w 4.4 x 10 12 /?6 (l + 0.024pg 2/3 T 6 ) erg K^cm" 3 . (B.13) 
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For estimates of the thermal conductivity we rely on the results obtained by Schaaf 
(1988, 1990). He investigated the cooling of a neutron star with magnetized envelopes 
and therefore calculated carefully the transport coefficients in the density range p < 10 10 
g cm -3 for a temperature of about 10 6 K and magnetic field strengths of 10 12 — 10 13 
G. The electron-ion collisions are insignificant as transport processes, because the iron 
ions form a crystal lattice in the considered temperature-density range. Since in the 
outermost layers of the neutron star the heat is transported both by radiation and by 
conduction, k = K ra d + K con d, we have to estimate first the relative importance of these 
two cooling mechanisms. 

For densities p > 10 3 g cm~ 3 , the radiative transport in the temperature and magnetic 
field range under consideration is limited by the opacity due to free-free transitions, 
k w kff (see Fig. 3 of Schaaf 1990). For p > p s {Bi2 = 1) the component of h/f parallel 
to the magnetic field is about 10 3 g _1 cm 2 . The corresponding radiative heat conductivity 
Krad = 16(tT 3 /(3k/9 s ) « 6.8 x 10 erg cm _1 s _1 K _1 . This value varies not much with 
increasing B12 because the growth of p s compensates at least partially the decrease 
of k. This K ra d has to be compared with the heat conductivity limited by electron- 
phonon collisions; the contribution of electron-impurity collisions is negligible for such 
low densities. For densities below 10 6 g cm -3 the parallel component of the conductivity 
tensor depends strongly on the local magnetic field strength because at low densities 
the quantization effects on the electron motion are dominant (Schaaf 1988, Fig 4.13); 
K CO nd(p s (B 12 = 1)) « 1.6 x 10 11 erg cm" 1 s" 1 K" 1 while K cond { Ps {B 12 = 10)) « 1.6 x 10 12 
erg cm -1 s _1 K _1 . Therefore, the heat transport in the polar cap region is nearly solely 
determined by its electronic part, k « K con d- 

Now we can estimate the e-folding cooling time in the magnetized outermost surface 
layer of the polar cap. Using B.1.8 for the depth of heat deposition and assuming that 
the heat by the bombardment of the polar cap surface is released in a depth of about 10 
radiation length (see CR80) we obtain 



Therefore, the strong local surface magnetic field is also the reason for an extremely short 
characteristic cooling time in the polar cap surface layer. For the thermostatic regulation 
described in this paper, the cooling by a few percent of 10 6 K may proceed in a time 
interval as short as 10 ns. Note that the independent estimate of the cooling time as 
given by CR80 is 



which for values of k and C given above, B 12 = 10 , T = 10 6 K and T - T = 10 4 K 
yields a cooling time t coo i w 20 ns, a few percent of our e-folding time r coo /. 



T CO oi(Bi2 = 1) « 200/is and t coo i(B 12 = 10) 



1/iS 



(B.14) 




(B.15) 
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B.3. Heating timescales 

The inflow of energy per unit time from the bombardment of the polar cap surface with 
ultrarelativistic electrons (positrons) ?77acc me c 3 uqj must be balanced by the change of 



Lorentz factor 7 acc is given by Eq. (A. 13) and depends on the shielding factor rj as well 
as on the local temperature, magnetic field strength and curvature. The characteristic 
heating time can be defined by 

C LT , 

Thect = 3 • (B.16) 

r) 7acc m c c J noj 

We will give a lower limit for Theat , i-e. consider the situation rj — 1 , when the accelerating 
potential drop is not yet shielded by thermionic ions (electrons). Using the expressions 
derived in Appendix A.l. and B.l. we find the characteristic heating time 



which returns for the B 12 = 10, p e = 0.07, L_ 3 = 2, Pi = 1, TZ 6 = 0.1 and T 6 = 1 a 
characteristic heating time Theat w 250 ns. Note that actually Theat oc r/~ 2 and for pulsars 
considered above r\ < 0.2. Therefore, the real characteristic heating time should be of the 
order of 10 /is, which is about (10 — A0)h/c, as determined by RS. 

Appendix C: Electron flow (ft • B > 0) 

Akin to the case f2 • B < 0, where positively charged iron ions are pulled out from the 
surface to contribute to the GJ density in the acceleration region, for the other "half" 
of pulsars, frequently called as "antipulsars" , the situation with £1 • B > has to be 
considered. The question is: what is the threshold surface temperature T s (in case of the 
thermionic particle extraction) or - electric field (in case of the field emission particle 
extraction), below which the binding energy of the electrons in the polar cap matter is 
significant to allow for the establishment of the charge depleted acceleration region. 

C.J. Thermionic electron emission 

In order to screen totally the vacuum gap electric field Eu the thermionic extracted 
electrons have to provide the Goldreich-Julian charge density. Since the electrons are 
nearly instantenously accelerated to relativistic energies, the current density due to this 
thermionic emission jth must somewhat exceed the Goldreich-Julian current density, i.e. 

3th > n-cje c, (C.l) 

(see Usov & Melrose 1995). The thermionic current density is described by the 
Richardson-Dushman equation (e.g. Gopal 1974) by 



the internal energy within the polar cap volume cx r: 



p L, which is CL dT/dt, where the 



r heat « 7.2 x 10- 8 p 6 T 6 5/3 L_ 3 PiX 4 ' 




s 



(B.17) 




(C.2) 
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The work function w has to exceed somewhat the electron Fermi energy ep fZiman ll972[) . 
Since it is a very complicated task to calculate w for the electrons at the neutron star 
surface, a task which has not yet been solved, we will use w ~ ep as a reliable lower 
limit. 

In the outer layers of the polar gap where the matter density is much lower than 
10 6 g cm~ 3 , and the temperature <C 10 7 K, the electrons are degenerated but non- 
relativistic, i.e. their Fermi energy ep is related to the Fermi-momentum Pp by 

p2 

£F = w^, (C3) 
2m*' y ' 

where for densities < 10 6 g cm -3 the effective electron mass to* coincides almost exactly 
with its vacuum value m e . Under conditions of the magnetized polar cap, where the 
magnetic field strength is of the order of 10 13 G, up to densities of ~ 10 5 g cm~ 3 only 
the ground Landau level (n = 0) is populated. In that case the Fermi-momentum of the 
magnetized electron gas is given by Fushiki et al. (1989) 

JM» = o) = ^p5. (C4) 

en 

Correspondingly, the electron Fermi-energy reads 

, 2TT 4 h i C 2 Z 2 p 2 . . 

C '(»=0)= e 2 mlmeA 2 B 2- (C5) 

Inserting p — p s (B) from Eq. (B.l), the electron Fermi energy at the solid iron surface 
of the polar cap is only a function of the magnetic field 

e F = 1.58 x lO- 9 B 2 i 5 erg « 986£^ 5 eV. (C.6) 

Approximating ep ~ w and with the Goldreich- Julian charge number density uqj = 
2ttB/ (P^TTce) = B/(P 1 ce), where Pi denotes the rotational period of the pulsar, we can 
rewrite Eq. (C.l) in the form 



2 . 78 x 10-1*^ < Tlexp - | } , (C.7) 



where T§ = T e /10 6 K and B\i = B s /10 12 G. This equation can be solved iteratively for 
given values of Pi and Bi 2 . By fitting dependencies on Pi and B12 we obtain an explicit 
expression for critical electron temperature (Eq. (CIO)) in the form 

T e ~ (4.5 x 10 5 )S?2 4 Pf 04 K. (C.8) 

For comparison see Eq. (10) in Usov & Melrose (1996) and Eqs. (2.13) in Usov & Melrose 
(1995), which we reproduced and refined here. The numerical factor 4.5 takes into account 
the most recent estimate of zero pressure surface matter density obtained by Lai ( 2001 ) 
for a 3 dimensional bec crystal (see Appendix B.l.). 
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C.2. Field electron emission 
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When the thermionic electron emission at the polar cap is not efficient, perhaps because 
of the fact that the neutron stars have been cooled down too much, still a field emission 
mechanism may provide a sufficiently strong electron flow. A non-zero En will cause a 
quantum mechanical tunneling of the electrons through the barrier provided by the work 
function and (e.g. Jessner et al. 12001(1 . The field emission current density of electrons is 
given by 

j f = ME {] exp{-N/E ]] }, (C.9) 

(e.g. Beskin 1982), where M = 3x 10 16 B 12 w-^ 2 s- 1 and JV = 2x 10 14 w 3 / 2 Vm _1 and w 
measures the work function in keV. Equating (as in the case of thermionic emission) the 
Goldreich- Julian outflowing current density with jf, and approximating again w ~ tp ~ 

2/5 

eV, the threshold electric field Ef,, above which the field emission of electrons 
becomes effective is given by 

1 < 3 x W 4 P 1 B^ /5 E ll {exp- (1.96 x 10 U B% 5 /E l{ )}. (CIO) 

This equation can be solved iteratively for given values of the pulsar period Pi and the 
magnetic field B 12 . For Pi = 0.1 and B X2 = 1 we find Ef,> 2.05 x lO^Vm^ 1 , a value 
certainly too large to occur at the polar cap surface (see also Appendix A in Gil & 
Melikidze 2002). 
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Table 1. Parameters of polar cap physics 
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0.17 


2.012 


2.025 


2.027 


3.9 


~ 10 


B0809+74 


0.17 


1.29 


11 


11 


> 150 
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~ 6 
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0.957 


0.94 


- 7 


B0826-34 
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1.85 


26 


~ 1 


~ 14 


14 


73 


33 


~ 8 


- 4.5 


~ 0.36 
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- 5 
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- 5 


0.071 


1.236 


1.239 


1.239 


7.9 


13.6 


B0031-07 


0.4 


0.94 


20 


- 6.8 


> 34 
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(1) P_i5 = P/10 -15 - period derivative 

(2) Pi - pulsar period in seconds 

(3) P2 - distance between driftbands in degrees of longitude 

(4) P3 = I//3 - primary drifting period in units of Pi (f^ - frequency of drifting feature) 

(5) P3 - circulational (tertiary) period in units of Pi (Eq. (A. 13)) 

(6) N = P3/P3 - number of circulating sparks (Eq. (A. 15)) 

(7) r p - polar cap radius in meters (Eq. (A. 4)) 

(8) d - circulation distance (s — d/r p < 1) 

(9) Vd — 2nd/ P3 - circulation speed in meters per second 

(10) r p /h - complexity parameter (h - gap height) 

(11) r, - shielding factor (Eqs. (5), (7) and (11)) 

(12) T s - actual surface temperature in 10 6 K (Eq. (10)) 

(13) T l - ion critical temperature in 10 6 K (Eq. (3) and (4)) 

(14) T e - electron critical temperature in 10 6 K (Eq. (5)) 

(15) B d - dipolar surface magnetic field in 10 12 G (Eq. (A. 3)) 

(16) B s - actual surface magnetic field in 10 12 G (Eq. (A. 2)) 



